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Abstract 

We derive a nonperturbative transverse momentum distribution 
for partons using a potential model to describe the quark-quark in- 
teraction inside the proton. We use this distribution to calculate the 
diferential cross-section of vr'^-production for intermediate values of 
transverse momentum in p — p collisions at high energies. Assuming 
a variable string tension constant for the quark-quark potential we 
obtain a very good description of the experimental data at different 
energies. The corresponding values of the mean transverse momen- 
tum of the partons are essentially lower than those obtained using a 
Gaussian transverse momentum parton distribution. 



1 Introduction 

During the last decades the tests of the perturbative QCD are focused on 
the experimental [1-15] and theoretical study [16-26] of hard processes like 
direct photon and 7r°-production with large transverse momentum in pp, pA 
and AA collisions. These processes offer a unique possibility to determine 
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the parton distribution functions (PDF) inside the proton as well as the par- 
ton fragmentation functions (PFF). In particular the good understanding 
of the pp data is the prerequisite for any attempt to extract new physics, 
related to the formation of a quark-gluon plasma phase, from the pA and 
AA data. Extensive studies of the vr (or 7) production in pp collisions have 
shown that the transverse momentum distribution glkx) of the partons in- 
side the proton has to be taken into account for a succesfull description of 
the observed pt spectrum [20-23]. In all these investigations one assumes a 
Gaussian form for gikx)- A new nonperturbative parameter is introduced 
through this approach: the mean instrinsic transverse momentum < kr > 
of the partons. Although the data of some experiments concerning hadron 
production at large pr [201 HI] could be explained with relative small < kx > 
values (~ 0.3 — 0.5 GeV), compatible with the Heisenberg uncertainty re- 
lation for partons inside the proton, there are a number of other processes 
leading to a large mean transverse momentum (< kx >~ 1 — 4 GeV), de- 
pending on Q^, for a description of the corresponding experimental data. 
Such a value of < /c^ > is too high and cannot be explained as an inter- 
nal structure of the proton j2S|- However, as mentioned by several authors 
1^ the form of gikr) can influence significantly the value of < fcy > as 
well as its pt dependence. In the present work we derive a transverse momen- 
tum distribution for the partons inside the proton using a potential quark 
model which has succesfuUy been used to describe the spectra of mesonic 
and baryonic |2H] bound states in the past. Following [2111211 investigate 
the three-body quantum mechanical bound state problem solving numeri- 
cally the Schrodinger equation and obtaining the single particle transverse 
momentum distribution for the consistuent parton. Our main assumption is 
that intrinsic transverse momentum effects are not influenced by the Lorenz 
boost along the beam axis and therefore could be treated within a nonrela- 
tivistic approach. We use the derived distribution to fit experimental data 
for the pp ^ TTo + X process. In particular we investigate the measurements 
for the pT-spectrum of the outcoming 7r° in three experiments performed at 
different center of mass energies. It turns out that a relatively low mean 
transverse momentum < /c^ > (pa 0(300 MeV)), compatible with intrinsic 
dynamics inside the proton, for the initial partons is sufficient in order to fit 
perfectly the experimental data. A smooth dependence of < /ct > on px, 
which within our approach is induced by a corresponding variation of the 
string tension in the quark-quark potential, is required. Our analysis shows 
that the fc^ distribution of the consistuent partons can be strongly infiuenced 
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by three-body effects and the form of the confining potential which have to 
be taken into account in order to describe correctly the experimental data 
concerning the pion production in pp collisions. The paper is organized as 
follows: in Section 2 we present the parton model differential cross section as 
well as the corresponding kinematics for the 7r°-production in pp-collisions. 
In section 3 we derive the intrinsic transverse momentum disrtibution for the 
partons inside the proton using the quark potential model of j2Z|- In section 
4 we present our numerical results concerning the description of the data of 
three different experiments 00113] as well as the corresponding dependence 
< kxipr) >■ Finally in section 5 we summarize our study and we discuss 
possible extensions of the present analysis. 



2 The pp ^ + X cross section 

In the lowest-order perturbative QCD (pQCD), the differential cross section 
for the direct pion production in pp collisions is given by: 

E^^ipp 7r° + X) = J dXadXbfa/p{Xa, Q'^)fb/p{Xb, Q"^) 



abed ' 



X —(ab^cd) — (1 



where fi/p {i = a, b) are the longitudinal parton distribution functions (PDF) 
for the colliding partons a and b while -Dtt/c is the parton fragmentation 
function (PFF) for the pion. For the scale Q we use the relation = 



2stu 



s2+P+u2 proposed in [12] with s, t, u the usual Mandelstam variables. The 
variable indicates the momentum fraction carried by the final hadron. 
The higher order corrections are taken into account by choosing K 2 
for the corresponding coefficient in (Q) in the pt region of interest. It is 
straightforward to include partonic transverse degrees of freedom using the 
following replacement [21] in the PDFs of eq.(|T|): 

dxi fi/p{xi,Q^) — > dxid'^kT,ig{kT,i) fi/p{xi, Q"^) (2) 

with i = a,b. In order to avoid singularities in the differential cross sections 
describing the partonic subprocesses we introduce a regularizing parton mass 
m = 0.8 GeV, as in ^H]; in the Mandelstam variables occuring in the de- 
nominator of the coresponding matrix elements. The explicit formulas of the 
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relevant partonic cross sections can be found in ^T]. Taking into account 
the transverse degrees of freedom we get the following expressions for the 
variables s, t, u: 

S = SXaXb H 2fcT,a " kr^b 

SXa^b 

t = -{Xa^ '-) \ kT,a-PT 

u = -(xfeH -) \ kT,b-PT (3) 

SXb Zc Zc 

Due to energy-momentum conservation the momentum fraction of the final 
hadron Zc is given by: 

Zc - . [V 

For a consistent description of the kinematics in the partonic subpro cesses 
we imply the cuts: 

Zc<l ; A;|_a < PtVs ; < pr^/s (5) 

To calculate the cross section given in eq.(^ we have first to determine the 
distribution gi^kx) and then perform the corresponding phase space integra- 
tions. Contrary to the usual treatment assuming a Gaussian form for g{kjn) 
we will here derive an alternative expression based on a widely applied quark 
potential model. 



3 The intrinsic transverse momentum distri- 
bution g(kT) 

Since the early days of quantum chromodynamics the main and almost unique 
tool used for the description of the hadronic bound states (mesons, baryons) 
remain potential models for the quark-quark and quark-antiquark pair inter- 
action. One of the most succesfuU models was proposed by A. Martin used 
first to describe mesonic states |2I] and later to describe baryons |2H]- The 
quark-quark interaction within this model is given as: 

V{r) = A,^r'-' + (6) 
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and a similar expression with adapted coefficients Agg, Bgg holds for the 
quark- ant iquark interaction. In the following we will consider a baryon con- 
sisting of 3 valence quarks interacting pairwise with the potential The 
Hamiltonian operator of the system is given as: 



2m 



(7) 



The baryonic ground state can then be obtained by solving the corresponding 
Schrodinger equation. Following ^28| it is convenient to introduce Jacobi 
coordinates 



6 
6 



r2 - ri 

2r3 - f2 - fi 

n + ^2 + ^3 



(8) 



getting the equation: 



2M 



€3 



m 



(v?, + v|j + v^(ei,6) 



(9) 

with M = 3m and m is the constituent quark mass. Eliminating the transla- 
tional mode (^3) we obtain a partial differential equation (PDE) depending 
solely on the variables ,^1, ^2- The potential energy V in eq.© is given as: 

= V{^,) + V{^ii^ - V36)) + n^(V36 +6)) (10) 



Introducing the radial variable ^ = yii+Q and the angular variables 
X, Oi,(pi,02, 02 we rewrite eq.(jn)) as follows: 



- -{^9^ (eO^^) + -\,-d^ (sin^ 2xdx^\ + ^ + - 



i sin 2x 



„ } 

cos^ X sm X 

V{i,xMAxMA2)^ = EG-^ (11) 



where the angular momentum operators are 

1 



9 (9 1 

sin^.M^^^"^ sin^^, 



^2 
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The reduced ground state wave function \^ can be expressed in terms of the 
hyperspherical harmonics Pl{^) as follows: 

^(e, X, ^1, 01, 02, 02) = E ^1, 01, O2, 02) (12) 

L=0 ^ 

The radial part of the wave function \& fullfils the equation |29j : 

15/4 + L(L + 4) m t> n n-?^ 

-Ml + T2^gUl - 72 tiL' v/ll' = (13) 



v 



As mentioned in [2H] the matrix Vl^/ = ^ dVLPl{Vt)V{i,VL)PL,{yL) for the 
particular choise of the potential © is diagonal dominant and therefore to 
a very good approximation the ground state of the baryonic system is deter- 
mined by the L = term in the expansion p2j) . Therefore the radial part of 
the ground state wavefunction obeys, within the above approximation, the 
equation: 

^ - + j2^Eo - V,,)u, = (14) 

where: Voo = ^oo + -Boo^°'^ with Aqq = f^gg and Bqq = \\Bqq. The constant 
A is given by: A = — r(1.55)r(1.5)r(3.05). It is straightforward to show 
that the ground state wavefunction in the momentum space (conjugate to 
the space of ^j, i = 1,2, 3) is given by: 

= N r die'^u,(i)^^Ml^^(M!!sO (15) 

J rC^ 

where J„ is the Bessel function of order n and iV is a normalization con- 
stant. It is useful to determine the transformation of the momenta k^^ to the 
cartesian momenta ki: 

4 = -^(^1 + ^2-24) 

% = -^ih + h + h) (16) 

The above expressions simplify in the center of mass frame k^.^ = of the 
baryonic system where the following relation is valid: 

kl = kl + kl + ki- k2 
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The two-particle density p{ki, ^2) in this case is given by: 



p{k„ h) = mki)\' = mki + ki + k, ■ k^)\^ (17) 

From eq.(|17p we obtain the one-particle transverse momentum density gik^) 



as: 



_^ roc rl poo _ y 

g(kT) = 47r dk, J ^dzj^ dk2kl\^{kl + kl + kl + zk2\Jkl + A;2)|2 (18) 

where z is the cosine of the angle between the vectors k = [kx, kz) and /c2- In 
fact the function Mo(0 can be obtained only numerically by solving equation 
fjl4p using the Numerov algorithm. Therefore also the transverse momentum 
distribution (fTHj) is known only numerically. The integrations in eqs. pSllSp 
can be performed to a great accuracy (relative error 10^^) using a mix- 
ture of Gauss-Kronrod quadrature and the VEGAS Monte-Carlo integration 
routine [201 • In Fig- 1 we present the density in transverse momentum space 
of a parton inside the proton obtained by our approach. The values of the 
constants Aqq and Bqq are choosen according to [2H1 in order to fit the size 
and the binding energy of the proton. The characteristic structure dominat- 
ing the form of g{kT) within our model is the second maximum at relative 
high transverse momenta. This leads, as we will see in the next section, to 
a reduction, relative to the Gaussian case, of the mean intrinsic transverse 
momentum of the partons needed to describe the experimental data for the 
7r°-production at various energies. 



4 Numerical results 

In the following we will use the derived distribution g{kT) in order to calculate 
within the parton model the differential cross section for the vr^-production 
in pp collisions according to eq.((H). For the longitudinal parton distribution 
functions we use the recent Martin, Roberts, Stirling and Thorne (MRST) 
scheme jSIj while for the parton fragmentation functions we use the Kniehl, 
Kramer, Potter [32] parametrization. The additional nonperturbative param- 
eter in our treatment is the mean value of the intrinsic transverse momentum 
< kx > which is related to the string tension Bqq in eq.®- Although the 
distribution g{kT) is derived for the valence quarks inside the proton here 
we will use the same distribution also for the initial gluons assuming uni- 
versality at the level of consistuent partons. In any case for longitudinal 



7 



momentum fraction Xi > 0.5 {i = a,b) the contribution of valence quarks 
is dominant and our description is accurate. The phase space integrations 
are performed using the VEGAS Monte-Carlo routine. In order to fit the 
experimental data we allow < kx > to vary as a function of the beam energy 
and the transverse momentum (px) of the finally produced hadron. We will 
analyse here the results of three experiments concerning 7r°-production at 
different energies. The first set of data are taken from the fixed target exper- 
iment performed in Fermilab (protons incident on H2 target) p. The cross 
section for the vr^-production with transverse momentum pt at midrapidity 
and for three different proton beam energies El = 200, 300 and 400 GeV 
is measured. In Fig. 2 we present the various datasets using symbols while 
with solid lines we show the results of our calculation and with dashed lines 
the corresponding results using a Gaussian glkx)- As we can see the two 
descriptions differ only in the region px < 1 GeV. Within our approach 
we need an intrinsic transverse momentum < kx > of the order of at most 
0.5 GeV in order to reproduce all the experimental data. It must be noted 
that in this case one can perfectly fit the data also for p^ < 1 GeV which, 
as we see in Fig. 2, is not possible using a Gaussian distribution gikx). In 
Fig. 3 we show the functions < kr > (pr) obtained using the quark model 
inspired function g{kT) as well as a Gaussian form. It is evident that the 
Gaussian model leads to much higher values oi < kx >. In particular this 
difference is large even for large values of px where our approach is more 
precise as the valence quark contribution is dominant. As already discussed 
in the previous section this difference relies on the fact that the non-Gaussian 
gikx) derived here posseses a second small local maximum at high transverse 
momenta (see Fig. 1) attributed to the form of the quark-quark potential and 
the many-body character of the system. 

The second experiment we consider is the WA70 at CERN SPS 6j. It 
is also a fixed target experiment with El = 280 GeV. We are interested 
in 7r°-production. In Fig. 4 we show the experimental data (full stars) and 
the corresponding pQCD calculation using the quark model inspired gi^kx) 
(solid line) as well as a Gaussian form (dashed line). Both distributions re- 
produce perfectly the experimental data and cannot be distinguished graph- 
ically. However our approach leads to significantly lower values of < kx > 
than in the Gaussian case. This can be clearly seen in Fig. 5 where we show 
the function < kxipx) > both for the quark potential model inspired g{kx) 
(full circles) as well as the Gaussian intrinsic transverse momentum distribu- 
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tion (full stars). Also in this case the difference in < > between the two 
approaches is large in a p^^-region where the main contribution is attributed 
to the valence quarks. For comparison we present also the same function for 
the Fermilab experiment at El = 300 Gev (open circles, see Fig. 3). 

Finally we have analysed the pp tt^ + X data of the most recent 
PHENIX experiment at the Relativistic Heavy Ion Collider with y/s = 
200 GeV pni- The corresponding cross section can be described to a good 
accuracy without inclusion of any /c^-smearing, a fact which is compatible 
with expectations for the perturbative character of the subprocesses involved 
in this case. Here we have fitted the experimental data using non-Gaussian 
/c^^-smearing effects. In this way we get a perfect description of the measured 
cross section. The obtained mean intrinsic transverse momentum is almost 
constant: < /c^ >~ 250 MeV. Based on Heisenberg uncertainty relation one 
could explain this value of the mean transverse momentum as an effect of the 
internal partonic structure of the incident proton. Our results are presented 
in Figs. 6,7. The PHENIX data (open cirlces) for together with the 

parton model calculations (crosses) using non-Gaussian g{kT) are shown in 
Fig. 6. The corresponding function < kx > (pr) is presented in Fig. 7. 

5 Concluding remarks 

Using a quark potential model capable to describe consistently baryonic sys- 
tems as three quark bound states we have derived an intrinsic transverse 
momentum distribution g{kT) of partons inside the proton. This, clearly 
non-Gaussian, distribution is characterized by the presence of a smooth local 
maximum at relatively high transverse momenta. Our approach is based on 
the idea that transverse momentum effects may be described nonrelativisti- 
caly as they are not influenced by the Lorenz boost along the beam axis. De- 
scribing the fcy-smearing effects in the parton model for pp collisions through 
this non-Gaussian distribution we calculated the differential cross section for 
7r°-production at midrapidity for different beam energies. Assuming that the 
corresponding nonperturbative parameter < kr >, related to g^kr), depends 
on pt of the finally produced hadron as well as the incident energy, we obtain 
a perfect description of the experimental data measured in pp collisions at 
three different experiments. The corresponding values of < fcr > as a func- 
tion of Pt could originate, according to Heisenberg uncertainty relation, from 
the internal partonic structure of the proton. Our approach is expected to be 
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valid also in single photon production as well as pA and AA processes ^3] . 
Therefore the performed analysis shows that many-body effects through a 
confining potential, reflected at the level of one-particle distributions, may 
influence strongly the fc-r-smearing phenomena observed in hadronic colli- 
sions and therefore should be taken into account for a correct description of 
the experimental data. 
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Figure 1: The intrinsic transverse momentum distribution 2n kTgikr) of a 
parton inside the proton obtained through the quark potential model de- 
scribed in section II. 
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Figure 2: The differential cross section for the vr^-pro duct ion in the Fermilab 
experiment The symbols represent the experimental datasets for three 
different beam energies. The solid lines represent the parton model results 
using a non-Gaussian 5'(/ct) while the dotted lines correspond to the analo- 
gous results using a Gaussian /cT-smearing. Only for < 1 GeV the two 
fits differ. 
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Figure 3: The functions < kripr) > for the three datasets presented in Fig. 2. 
We distinguish between the results of the parton model using non-Gaussian 
(n. G.) or Gaussian (G.) intrinsic transverse momentum distribution. 
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Figure 4: The differential cross section for vr^-production in WA70 at El = 
280 GeV (full stars) and the corresponding parton model calculation using 

— * 

a non-Gaussian (solid hne) as well as a Gaussian (dotted line) gikr)- 
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Figure 5: The dependence kripT) for the WA70 data set using a non- 
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Gaussian as well as a Gaussian g{kT) in comparison with the correspond- 
ing behaviour found for the Fermilab experiment at — 300 GeV with a 
non-Gaussian g{kT)- 
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Figure 6: The — > 7r° + X differential cross section for the PHENIX experi- 
ment {^/s — 200 GeV). The experimental data are presented by open circles 
while the crosses indicate the parton model calculation using a non-Gaussian 
/cT-smearing. Only the experimental errors can be seen at this scale. 
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Figure 7: The function < kxipr) > obtained using the non-Gaussian g{kT) 
in the parton model description of the PHENIX data. The sohd hne at 
250 MeV is shown as a guide. 
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